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This  report  is  one  of  a series  that  deals  with  the  problem  of  spinning 

flexural  dynamic  systems.  In  this  initial  document,  a general  overall 

description  of  the  hybrid  coordinate  approach  is  provided  as  an  introduction 

only.  The  general  structure  and  explicit  development  of  the  equations  have 
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been  provided  in  the  literature.  A hybrid  coordinate  approach  is  illus- 
trated by  applying  it  to  a simplified  spacecraft  system  so  that  the  general 
developments  appearing  in  the  literature  caii  be  better  appreciated  by  poten- 
tial users.  A simple  example  is  selected  to  illustrate  many  features  of 
flexural  systems  engaged  in  spin.  The  origins  of  specific  terms  in  the 
dynamical  equations  and  their  anticipated  effects  upon  the  dynamics  are 
discussed. 

A stability  analysis  is  provided  for  a simplified  spacecraft  system. 

The  emphasis  is  on  the  development  of  literal  stability  criteria  for  the 
wobble  motions  of  the  simplified  spacecraft  system  that  can  be  of  value 
in  preliminary  design  phases.  Admittedly,  some  risk  is  associated  with 
extrapolation  of  results  for  simplified  dynamical  systems  to  more  com- 
plicated excimples.  The  stability  criteria  can  be  used  to  duplicate  results 
of  previously  published  studies  involving  considerably  more  sophisticated 
models  as  a point  of  departure,  however.  The  stability  criteria  can  be  of 
value  in  preliminary  design  phases  but  are  not  intended  to  supplant  more 
detailed  followup  studies  to  verify  ultimate  system  stability. 

^P.  W.  Likins,  E)ynamics  and  Control  of  Flexible  Space  Vehicles.  Jet 
• Propulsion  Laboratory,  Technical  Report  No.  5Z-13Z9,  Revision  1, 

January  1970. 

2 

P.  W.  Likins,  Finite  Element  Appendage  Equations  for  Hybrid  Coordinate 
Dynamic  Analysis.  Jet  Propulsion  Laboratory,  Technical  Report  No. 
32-15Z5,  October  1971. 

3 

P.  W.  Likins  and  P.  H.  Wirsching,  "Use  of  Synthetic  Modes  in  Hybrid 
Coordinate  Dynamics  Analysis,  " AIAA  J.  , 6,  10:  1867-1872,  October  1968. 


A. 


HYBRID  COORDINATE  APPROACH 


Traditional  modal  analyses  entailing  modal  deformation  coordinates 
of  the  entire  vehicle  have  been  employed  for  certain  classes  of  spacecraft 
with  much  success.  These  approaches  possess  the  advantage  of  being 
rigorous  as  long  as  small  deformations  are  observed  and  can  be  computa- 
tionally efficient  in  that  irrelevant  high  frequency  response  can  be  truncated 
from  the  simulation.  Traditional  modal  analyses  entailing  modal  deforma- 
tion coordinates  of  the  entire  vehicle  may  not  be  applicable  to  a larger  class 
of  space  systems  frequently  encountered:  systems  of  time-varying  con- 
figuration, systems  employing  nonlinear  control  systems,  systems  employ- 
ing discrete  dampers,  and  systems  requiring  accommodation  of  large 
arbitrary  vehicle  motions. 

Faced  w"  se  obstacles,  some  analysts  have  resorted  to  discrete 

parameter  f'  s.  The  totcil  system,  including  flexural  appendages, 

is  idealized  collection  of  rigid  subbodies  articulated  to  each  other  in 

some  topological  tree  configuration.  Equations  of  motion  are  then  derived 
for  this  system.  In  this  approach,  such  obstacles  as  large  configuration 
variations,  inclusion  of  discrete  dampers,  nonlinear  controls,  or  large 
arbitrary  vehicle  motions  can  be  accommodated.  The  physical  system 
analyzed  may  not  be  particularly  representative  of  the  actual  system, 
however.  Satisfactory  simulation  may  require  many  rigid  bodies.  The 
analyst  in  this  case  is  faced  with  a complex  task  of  formulating  equations 
of  motion  for  many  rigid  bodies  and  is  forced  to  carry  along  irrelevant  high 
frequency  response  in  any  simulation  attempts. 

The  hybrid  coordinate  approach  represents  a natural  compromise  be- 
tween fully  discrete  and  exclusively  modal  methods.  The  most  general  defini-  * 
tion  of  the  hybrid  coordinate  approach  to  formulating  equations  of  motion  is 

any  approach  that  combines  discrete  coordinates 
describing  the  translations  and  rotations  of  some 
bodies  or  reference  frames  of  the  system  with  dis- 
tributed or  modal  coordinates  describing  the  small 
relative  motions  of  other  parts  of  the  system  .... 
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The  hybrid  coordinate  approach  is  applicable  when  all  or  part  of 
the  vehicle  admits  the  assumption  of  small  linearly  elastic  deformations. 
The  most  efficient  simulation  combines  discrete  coordinates  with  modal 
coordinates,  retaining  the  generality  of  discrete  coordinates  where 
necessary  and  securing  the  computational  advantages  of  modal  coordinates 
where  possible. 


The  literature  devoted  to  the  hybrid  coordinate  approach  received  a i 

major  impetus  in  197  0 with  the  publication  of  a JPL  report  entitled, 

Dynamics  and  Control  of  Flexible  Space  Vehicles,  by  P.  W.  Likins.^  At 
that  stage  of  development,  the  method  presented  was  formally  restricted 
to  systems  in  which  each  flexible  appendage  is  attached  only  to  a rigid 
body  or  to  several  rigid  bodies  interconnected  so  that  relative  motion  can- 
not induce  deformation  of  the  flexible  appendage.  The  method  as  presented 
is  best  suited  to  systems  in  which  the  flexural  appendage  is  the  outer  member 
of  a system  of  interconnected  bodies  of  topological  tree  configuration.  Any 
number  of  flexural  subsystems  could  be  accommodated  subject  to  the  restric- 
tions, and  furthermore  the  appendage  could  be  engaged  in  spin.  j 


The  dynamics  of  these  general  spacecraft  configurations  were  derived 
using  an  Eulerian  formulation.  Generalized  matrix  expressions  were  devel- 
oped and  presented.  These  formidable  general  expressions  can  be  specialized 
to  particular  cases  as  long  as  they  fall  within  the  stated  restrictions.  Since 

the  publication  of  Ref.  1,  several  contributions  have  developed  the  idea  of 
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a hybrid  coordinate  approach  even  further. 


B.  SCOPE  OF  WORK 


The  following  effort  parallels  to  some  extent  the  work  of  Likins  and 
5 

Barbara.  It  is  hoped  that  this  report  might  serve  as  a first  step  toward  a 


4 

P.  W.  Likins  and  G.  E.  Fleischer,  Large -Deformation  Modal  Coordinates 
for  Nonrigid  Vehicle  Dynamics,  Jet  Propulsion  Laboratory,  Technical 
Report  No.  32-15^5,  November  1972. 

5 

P.  W.  Likins  and  F.  J.  Barbara,  Attitude  Stabili^  of  Spinning  Flexible 
Spacecraft,  UCLA  Engineering  Report  No.  7176,  Decenri^ber  1971;  also 
NASA-CR-123569. 
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comprehensive  understanding  of  the  dynamics  of  spinning  flexural  systems. 

In  this  document,  essential  features  of  flexural  systems  engaged  in  spin  are 
discussed.  Since  we  aim  to  develop  a general  awareness  of  the  more  impor- 
tant features  of  flexural  dynamical  systems,  the  discussion  is  directed  toward 
the  nonexpert  in  this  area.  Hopefully,  this  report  and  subsequent  documents 
will  serve  as  a palatable  transition  to  the  literature  in  the  field. 
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II.  DYNAMICAL  FORMULATION 


In  this  section,  a particularly  simple  example  is  analyzed.  This 
example,  despite  its  simplicity,  reveals  many  features  of  flexural  systems 
engaged  in  spin.  The  origins  of  specific  terms  in  the  dynamical  equations 
and  their  anticipated  effects  upon  the  dynamics  are  discussed. 

The  example  consists  of  a rigid  hub  and  a pair  of  mass  particle/ 
massless  elastic  spring  combinations  mounted  symmetrically  on  the  hub. 

The  mass  particle/massless  elastic  spring  combinations  are  assumed  to 
be  representative  of  the  flexural  subsystem.  The  nature  of  th  intercon- 
nection between  the  flexural  subsystem  and  the  rigid  portion  is  optional. 

A variety  of  different  interconnections  are  assumed;  radially  mounted  sub- 
systems, anti  cantilever  subsystems,  and  an  intermediate  combination  of  the 
two  extremes.  We  will  show  that  the  nature  of  the  interconnection  can  be 
important  and  can  influence  the  dynamics  in  a predictable  way.  It  is  further 
assumed  that  this  entire  system  is  engaged  in  spin.  The  nominal  spin  state 
is  pure  spin  about  a symmetry  axis  of  the  vehicle.  This  system  is  illustrated 
in  Figure  1. 

The  equations  of  motion  for  the  appendage  are  derived  in  Section  II.  A. 
The  effects  of  spin  are  discussed  in  Section  II.  B.  In  Section  II.  C,  the 
equations  of  motion  for  the  entire  system  are  derived,  based  on  the  radial 
mounting  of  the  flexural  subsystem  to  the  rigid  hub.  Again,  important 
properties  of  the  overall  system  are  discussed. 

A.  DYNAMICAL  FORMULATION  OF  APPENDAGE  SUBSYSTEM 

The  simple  model,  shown  in  Figure  2,  consists  of  a rigid  core  to 
which  two  particles  are  attached  through  springs.  The  principal  axes  of 
the  system  in  its  nominal  state  remain  coincident  with  the  principal  axes 
of  the  core,  as  the  particles,  each  of  mass  m,  are  symmetrically  located 
a distance  | r.]  along  the  axis  y.  The  particles  are  allowed  to  move  radially 
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along  the  y axis,  as  well  as  rotationally  through  a 2-degree-of-freedom 
pivot  at  the  attachment  point,  thus  allowing  deformations  in  the  other  two 
axes  as  well.  Conceptually,  we  might  visualize  the  particle  as  being 
constrained  to  move  along  the  massless  tube.  Motions  of  the  particle 
transverse  to  the  tube  direction  are  achieved  by  a 2-degree-of-freedom 
hinge  located  at  the  interconnection  with  the  spacecraft.  The  notation  to 
be  used  in  the  following  text  is  indicated  in  Figure  Z, 

Let  jr  j and  2 the  equilibrium  orientations  of  the  respective  mass 
particles  when  the  system  is  engaged  in  a state  of  pure  spin  about  the  sym- 
metry axis.  Let  and  62  represent  small  local  deformations  of  the 
respective  particles  from  their  equilibrium  positions.  The  translational 
equations  of  motion  for  the  general  particle  i can  be  obtained  by  application 
of  Newton's  second  law  of  motion,  which  states 


N,2 

m. (X.)  = E f' 
dt  J J 


(1) 


Nd2 

where — r-(  ) denotes  a second  inertial  time  derivative  of  the  quantity  in 
dt 


parentheses,  and  j is  an  arbitrary  siunmation  index, 
contact  forces  acting  on  particle  i is  given  as  E _F^. 

j ^ 


vector  X.  into  component  parts,  we  have 


The  sum  of  all 
Breaking  down  the 


N,2 

m.  —5-  (R  + D + r.  + 6.) 
‘ dt^  - ^ 


(2) 


By  definition  of  the  mass  center,  we  have 


mT£^  +Em.  (r. 
1 


+ ii) 


= 0 


(3) 
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By  symmetry  and  equality  of  masses,  we  have 


Em.£.  = 0 


Thus,  the  expression  for  reduces  to 


E m.  6 . 

1 — i 


(4) 


Further  assuming  that  we  are  dealing  with  a force-free  system  yields 


-^R  = 0 (5) 

dt 


As  a result  of  these  assumptions,  Eq.  (2)  can  be  expressed  as 


N,2 

m.  — =-  (6. 
■ dt^ 


2 m. 


+ 


ii) 


(6) 


where  ? Fj  represents  the  vector  sum  of  all  contact  forces  acting  on  par- 
ticle i.  Again,  it  is  assumed  that  there  are  no  net  external  forces  acting 
upon  the  total  system. 


A A A T 

Define  an  orthonormal  triad  of  unit  vectors  (x,  y,z)  that  are  fixed 
in  the  rigid  hub.  The  inertial  angular  velocity  of  this  basis  is  designated 
by  u)^.  The  vectors  r.  are  invariant  with  respect  to  this  basis,  but  the 
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vectors  are  not  invariant  in  this  basis.  Carrying  out  the  indicated 

i 5!< 

differentiation  by  successive  application  of  the  Coriolis  theorem,  we  obtain 


) = )+„x(  I 


(7) 


and  letting  m = = ^2*  find  that  Eq.  (6)  reduces  to 


( 


00 


m ;6.v2wx6.+d)x6.+wx(wx6. ) 
I -1  — -1  - -1  - - 1 


£.  I-  00  U -I 

- L -22-  6.  + 2wx6.+wx6.+cjx{wx6.) 

j_  j L J J j J J 


wxr.+wx(uxr.)(  = S F / 
- -X  - -X  j j -J 


(8) 


where  (°)  denotes  differentiation  with  respec^to  a rotating  frame  of  reference 
and  is  a shorthand  notation  for  the  operator  ( ).  In  the  following,  we  are 


interested  in  developing  perturbation  expressions  about  a nominal  case  of  pure 
spin.  By  application  of  the  angular  velocity  chain  rule,  we  can  write 


0)  = S2  + 6 u 


(9) 


where  is  the  nominal  angular  velocity  of  the  system  in  its  equilibrium 
state  of  pure  constant  spin.  The  terms  retained  in  the  following  expressions 
are  consistent  with  first  approximation  analysis.  Products  of  small  quantities. 


The  operator  -gp-  ( ) designates  differentiation  of  the  quantity  in  parentheses 
with  respect  to  a rotating  vector  basis  P.  The  inertial  angular  velocity  of 
that  basis  is  a>. 
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PBB^ 
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such  as  deformations  and  variational  angular  velocities,  are  considered 
second-order  small. 


The  first-order  approximation  to  Eq.  (8)  is 


+ 2fi  x^.  - 2 L — fix 

j “‘>j'  j ^ j J 


m 


[oo  ^ oo 

Li  - 2 -^li 

1 j m.j,  j 

m.  1 . 

+ n X ( S2  X 6 . ) - 2-;^  S2  X ( S2  X 6 . ) - S' 

— ' — — 1 j m.p  — “ “J  J ^ 

^ i + Lli  ^ (0.  * Hj) 

+ nx  (6oj  X r +F'^  _i  = 1,  2 


-m 


:io) 


where  the  readily  identifiable  homogeneous  terms  in  deformations  have 
been  grouped,  and  svimmations  over  j are  over  j = 1 and  2. 

A portion  of  the  contact  forces  attributable  to  the  elastic  restoring 
forces  has  been  partitioned  and  grouped  on  the  left-hand  side.  This  con- 
tribution is  considered  proportional  to  the  deformations  from  the  equilib- 
rium configuration,  and  the  "weighting"  matrix  that  maps  these  deformations 


into  an  elastic  force  F is  referred  to  as  the  stiffness  matrix.  The  remain- 


—  e 


ing  contribution  to  the  contact  forces  is  designated  F'^  and  is  inferred  to  as 
preload  effects. 


The  equation  formulation  is  encumbered  by  the  presence  of  terms 


that  reflect  the  effects  of  mass  center  shifts,  i.e.,  )•  ® 


simple  transformation  of  variables,  those  terms  involving  the  mass  ratio 
can  be  removed  in  a rigorous  treatment.  In  fact,  such  an  approach  has 
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much  value  because  it  emphasizes  some  important  features  about 
flexural  systems.  Consider  the  transformation 

i.=Ai  -±2 

JL  =Ai  +±2 

Applying  this  transformation  pair  to  Eq.  (10)  yields  Eqs.  (12)  and  (13). 

By  differencing  the  two  vector  Eqs.  (10)  and  making  the  identification 
in  Eq.  (11),  we  have 

TOO  o I 1 2 

m [ + 2n  X + 2 X (n  X |)J  “ 

= -2m  ^ ^ (II  ^ + 6a>  x (n  x ^^) 

+ n X (^x  r^)]  + F'^  - F'^  (12) 

and  by  adding  the  two  vector  Eqs.  (10)  and  making  the  identification 
in  Eq.  (11),  we  have 

[oo  oo  o o 

1 - 2^  TL  + 2£  x 7]_  - 4^  n x J 

+ 

= F'^  + F'^  (13) 

Two  observations  are  made.  First,  the  terms  attributed  to 
mass  center  shifts  are  absent  in  Eq.  (12).  Second,  only  Eq.  (12) 


n X (n  X 11)  - 2 -22-  nxinxii)  - f^-f'^ 

— — m.p  — — -L'  — e — « 
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contains  forcing  terms  that  reflect  general  spacecraft  motion.  If  we 
further  identify  the  correspondence  of  the  variable  to  an  asymmetric 
mode  and  the  variable  to  a symmetric  mode,  then  the  observations 
take  on  new  meaning.  Symmetric  mode  shapes  are  not  likely  to  be 
excited  by  rotational  motions  of  the  hub.  Conversely,  an  excited  sym- 
metric mode  does  not  usually  produce  rotational  motions  of  the  hub. 

The  absence  of  coupling  with  spacecraft  rotational  motions  is  consistent 
in  this  example.  Another  appealing  fact  is  that  terms  attributable  to  mass 
center  shifts  appear  only  in  Eq.  (13).  Mass  center  shifts  do  not  occur 
with  asymmetric  mode  shapes.  Section  II.  B.  discusses  the  effects  that 
arise  primarily  as  a consequence  of  spin.  Equation  (12)  in  particular  is 
dealt  with  at  some  length. 


B.  EFFECTS  OF  SPIN 

The  left-hand  portion  of  Eq.  (12)  is  repeated  for  convenience: 


Coriolis  Centripetal 


(14) 


Two  terms  in  particular  arise  as  a consequence  of  the  dynamics 

formulation:  Coriolis  and  centripetal  acceleration  effects.  To  better 

appreciate  the  role  of  these  effects,  we  coordinatize  the  equations  in 

^ ^ ^ T 

a basis  fixed  in  the  hub  (x,y,z)  . 

The  nominal  inertial  angular  velocity  of  the  system  in  its  equilibrium 
orientation  is  given  by 


n = n z 


(15) 


We  anticipate  that  the  preload  effect  is  not  functionally  related  to  hub 
motions  in  Eq.  (13);  see  Section  n.B.2. 
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Equation  (14)  can  be  rewritten  as  follows.  The  vector  cross- 
product  operation  can  be  expressed  as 


axb  = {x}^ab  = 


^2' 

-aj 

0 


(16) 


Using  this  isomorphism,  we  can  write  the  vector  expression  as 


{x}'^[m  {1}  + 2m[n]  {i}  +m  [n][nl{|}  + [K]{|}j  (17) 

where  the  difference  in  elastic  forces  has  been  replaced  by  [k]  {!}• 


1.  CENTRIPETAL  ACCELERATION  EFFECTS 

The  contribution  from  the  centripetal  acceleration  term  is 


0 o' 

m[n][n]{|}  =m 

0 0 

. 0 0 0. 

§ 

(18) 


The  deformation  dependent  terms  collectively  become 


fk  -mn  0 01 

X 

0 k -mQ^  0 

§ 

y 

0 Ok. 

z 
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where  for  convenience  the  stiffness  matrix  was  assumed  diagonal.  The 
general  effect  of  the  centripetal  acceleration  is  a reduction  in  stiffness 
or,  analogously,  a reduction  in  the  unloaded  frequency  of  oscillation  of 
the  structure.  "Unloaded"  denotes  that  the  effects  of  preload,  a 
deformation-dependent  effect,  have  not  been  incorporated.  The  pre- 
load effect  is  still  retained  in  the  contact  forces  F'^.  Preload  effects 
and  the  elastic  restoring  effects  collectively  result  in  a stiffness  matrix 
that,  when  suitably  normalized,  yields  a set  of  loaded  frequencies. 

These  loaded  frequencies  are  altered  by  the  centripetal  accelerations 
as  discussed  previously. 

Terms  in  deformation  rates  also  appear.  Again,  this  is  a spin- 
induced  effect.  The  structure  of  the  equations  is  altered  by  the  presence 
of  these  terms  from  the  form  involving  only  mass  and  stiffness  matrices 
to  the  more  general  homogeneous  form 

Mq  + Gq  + Kq  = 0 

The  gyroscopic  contribution  Gq  has  a unique  structure  in  that  the  G 
matrix  is  skew  symmetric.  The  role  of  the  Coriolis  acceleration  can 
be  important,  but  this  discussion  will  be  deferred  until  the  remaining 
contact  forces,  called  preload  effects,  are  evaluated. 

2.  PRELOAD  EFFECT 

The  remaining  forces  are  evaluated  for  three  special  cases:  a 
radially  mounted  subsystem,  an  anticantilever  subsystem,  and  an 
intermediate  combination  of  the  two  extremes.  These  configurations 
are  illustrated  in  Figure  3.  Anticipating  the  results,  we  find  some 
marked  differences  between  the  configurations.  The  motivation  for 
considering  these  three  excimples  was  provided  in  Ref.  5.  The  re- 
maining contact  forces  in  the  radially  mounted  case  are  determined 
as  follows.  In  its  nominal  equilibrium  configuration,  the  massless 
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spring  supporting  the  mass  particle  is  in  tension,  the  magnitude  of 
2 

which  is  mn  r , where  r could  be  further  broken  down  into  a nominal 

y y 

unloaded  position  plus  a contribution  induced  by  spin.  For  convenience 
of  exposition  only,  assume  that  the  mass  flexural  subsystem  attach- 
ment point  to  the  rigid  hub  is  at  the  point  0.  This  configuration  is 
illustrated  in  Figure  4. 

To  first  approximation  analysis,  the  remaining  contact  force  F'^, 
resolved  into  the  coordinate  basis  [x,  y,  z]  , is 


- m n 6 


F'^  = {x}^  < - m r 


(21) 


- m n 6 


where  the  contribution  from  the  preload  effect  tends  to  provide  a 
restoring  force  and  as  such  would  positively  augment  the  stiffness 
matrix  [ deformation- dependent  terms  in  Eq.  (12)].  A general 
stiffening  effect  is  achieved  through  the  preload  effect  for  the  radially 
mounted  configuration.  This  effect  is  configuration  dependent,  as 
shown  in  the  following  example. 

Consider  the  case  of  the  anticantilever  particle  illustrated  in 
Figure  5.  Again,  assume  that  the  particle  in  its  nominal  equilibrium 
position  is  displaced  to  a distance  r^,.  For  convenience  it  was  assumed 
that  the  base  of  the  spring  is  connected  to  the  satellite  at  point  P,  a 
distance  of  2r  from  the  nominal  center  of  mass.  The  interconnect- 

y 

ing  spring  is  in  compression.  A better  pictorial  representation  of 
the  nature  of  the  anticantilever  mount  is  presented  in  Figure  3. 


The  application  of  a transformation  [Eq,  (11)]  does  not  alter  this 
basic  conclusion. 
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The  contact  force,  when  resolved  into  the  coordinate  basis  [x,  y,  z]  , is 
to  first  approximation 


F'^  = [x]"^ 


m 6 


- m fJ  r 


2 

m n 6 


(22) 


z -■ 
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The  components  of  force  that  are  proportional  to  deformations  are  not 
restoring  in  this  case.  Thus,  the  stiffness  matrix  is  augmented  such 
that  a general  reduction  in  stiffness  can  be  observed  due  to  preload 
effects.  We  can  conclude  that  preload  effects  are  configuration  dependent. 
One  final  example  reillustrates  this  fact. 

Consider  the  intermediate  case  where  the  particle  is  neither  radially 
mounted  or  anti  cantilever.  The  only  preload  effect  in  this  instance  is  the 
steady- state  force  required  to  maintain  the  equilibrium  orientation  of  the 
particle  and  is  directed  along  -y.  In  this  case,  preload  effects  do  not 
contribute  to  the  stiffness  matrix,  whereas  they  altered  the  stiffness 
matrix  in  the  previous  cases.  In  the  radially  mounted  case,  there  was 
a general  stiffening  due  to  the  preload  effects.  In  the  anticantilever  case, 
there  was  a general  reduction  in  stiffness  due  to  the  preload  effects.  These 
facts  are  summarized  in  Table  1. 


Table  1.  Effect  of  Preload  on  Stiffness 

Properties  for  Varying  Configurations 


Configuration 

Qualitative  Effect 
on  Stiffness 

Quantitative  Effect  on  Stiffness^ 
/Loaded  Frequency\ 

\ Squared  / 

Radially  Mounted 

Increased  stiffness 

Eq.  (21 );  k/m  + Cl^ 

Anticantile  ve  r 

Reduced  stiffness 

Eq.  (22);  k/m  - 

Neutral 

No  effect 

a 1/2 

Unloaded  frequency  (k/m)  , when  augmented  by  preload  effect,  is  defined 

as  loaded  frequency. 
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Collectively,  the  effects  of  preload,  nominal  elastic  restoring  force 
effects,  and  centripetal  accelerations  are  summarized  in  Table  2. 


Table  2.  Sximmary  of  Effects  of  Preload  a^d  Centripetal 
Acceleration  on  Stiffness  Matrix 


/ Loaded 
yFrequency/ 

Radial 

Anticantilever 

Neutral 

2 

a 

X 

^x  2 2 

(-J  + n^)  - 

m 

(—  - n^)  - 

m 

k y 

- n 

m 

2 

a 

y 

k - 

{-^)  - 
m 

k y 

- n 

m 

k y 

° 

2 

a 

z 

k y 

(-^  + n ) 

m 

k , 

(—  - n^) 

m 

k 

Loaded  frequencies  resvilting  from  collective  effects  of  elastic  restoring 
forces  and  preload  effects  are  in  parentheses.  Contributions  from 
centripetal  acceleration  effects  appear  outside  the  parentheses. 

3.  CORIOLIS  EFFECTS 

We  pointed  out  that  terms  in  deformation  rates  appeared  in  Eq.  (14) 
and  that  they  were  clearly  a spin- induced  effect.  The  basic  structure  of 
the  equations  is  altered  by  the  presence  of  these  terms  from  the  form 
involving  only  mass  and  stiffness  matrices  to  the  more  general  homogeneous 
fo  rm 

Mq  + Gq  + Kq  = 0 (23) 

Further,  the  gyroscopic  contribution  Gq  has  a unique  structure  in  that  the 
G matrix  is  skew  symmetric.  The  role  played  by  this  effect  is  examined 
in  the  following  paragraphs. 
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The  homogeneous  portion  of  Eq.  (14)  specialized  for  the  case  of 
radially  mounted  particles  is 

[M|  + 2M0  i + K'  |]  = 0 


where  the  matrix  K'  reflects  the  contributions  from  centripetal  accelera- 
tion effects,  preload,  and  normal  elastic  restoring  forces.  The  scalar 
equivalent  to  Eq.  (24)  after  dividing  through  by  the  mass  is 


Transverse 
to  Spin 


§ -2n  § 

X y 


+—  § = 0 
m =x 


+-^  1=0 
m ^y 


I 


+—  § = 0 
m 


Coriolis  Centripetal  Preload  Elastic 

The  equations  of  deformation  decouple  into  equations  transverse  to  the 
spin  direction  and  along  the  spin  direction.  Consider  the  equations  of 
motion  transverse  to  the  spin  direction  (out-of-plane) 


where 


L - 2n  § + a 5 =0 

^ y X ^x 


§ +201  +a^§  =0 

^y  X y ^y 


ok  o " 

.,2  X j 2 y 

c = and  a = — i-  - O 

X m y m 


If  the  hub  were  vanishingly  small,  then  the  frequencies  considered  here 
would  be  system  frequencies.  If  we  include  the  effects  of  the  hub,  the 
result  of  the  Coriolis  coupling  effect  might  be  altered.  With  this  restric- 
tion, if  we  transform  Eq.  (26)  into  the  S domain  and  evaluate  the  charac- 
teristic equation,  we  have 

(S^  + (S^  + a^)  + = 0 

X y , - _ - 

Coriolis 

Coupling 


(S^+  <J  )(S^  + a‘) 
' X y 


In  Figure  6 the  poles  represent  the  locations  of  frequencies 

a and  a . They  are  located  on  the  complex  axis  since  we  assume 
X y 

there  is  no  damping  in  the  system.  The  Coriolis  coupling  effects 

are  represented  by  the  factor  K.  The  inclusion  of  Coriolis  coupling 

can  alter  the  system  frequencies  but  does  not  alter  the  fact  that  the 

''system"  roots,  irrespective  of  gyroscopic  coupling,  must  be  on 

the  imaginary  axis  in  the  event  of  no  damping.  The  net  effect  is  a 

further  separation  of  the  loaded  frequencies.  The  highest  root  is 

pulled  still  higher,  while  the  lowest  root  is  further  pulled  to  a lower 

frequency.  Our  intuition  can  be  important  in  assessing  the  impact 

of  this  coupling  effect.  In  the  case  of  a spinning  helicopter  blade, 

we  might  assume  the  blade  would  be  nearly  inextensible  radially. 

The  loaded  frequencies  of  the  substructure  would,  in  this  case,  be 

2 2 

essentially  unaltered  by  the  Coriolis  coupling  effect  ). 

In  general,  the  importance  of  the  Coriolis  coupling  effect  must  be 
assessed  in  each  physical  application  and  should  not  be  casually 
dismissed. 
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c. 


COMP OSITE  EQUATIONS  OF  MOTION 


We  have  concentrated  primarily  on  the  development  of  equations 
of  motion  for  the  Qexural  subsystem.  To  complete  the  system  descrip- 
tion, we  must  augment  this  subset  with  equations  of  motion  for  the 
entire  spacecraft. 

This  general  approach  has  been  delineated  in  Ref,  1 for  more 
general  dynamical  systems  than  those  presently  under  consideration. 

In  particular,  the  general  structure  of  equations  is  developed  for  the 
following ; 

1.  Flexural  appendages  attached  directly  to  a rigid  hub 

2.  Flexural  appendages  engaged  in  spin 

3.  Systems  incorporating  discrete  dampers  and  control 
system  nonlinearities 

The  general  structure  of  the  equations  developed  in  Ref,  1 could 
be  used  as  a starting  point  for  the  simple  example  described  here.  The 
approach  adopted  here  involves  a development  of  composite  equations 
of  motion  from  first  principles  and  parallels  to  some  extent  the  develop 
ment  of  equations  in  Ref.  1.  However,  the  equation  development  is 
specialized  to  the  example.  Hopefully,  this  might  serve  as  an  intro- 
duction to  the  more  comprehensive  development. 

In  the  following  paragraphs,  the  composite  equations  of  motion 
are  developed.  Since  our  aim  is  the  development  of  expressions 
suitable  for  linearized  stability  analysis,  a variational  approach  is 
adopted  and  first  order  approximation  expressions  are  derived. 

The  total  angular  momentum  of  the  composite  system  referenced 
to  the  composite  center  of  mass  is  defined  as 

H*^  = r £ X £ dm 
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where  the  vector  £ emanates  from  the  instantaneous  center  of  mass 
and  terminates  at  a differential  mass  particle  dm.  The  operator  (*) 
is  understood  as  an  inertial  time  derivative,  whereas  the  operator  (°) 

denotes  a time  derivative  with  respect  to  a designated  reference  frame. 

c 6 

11  can  alternately  be  expressed  as 


h'"  = °J  . w + m^  ^ X £^  f-y*  r X r 


dm 


(30) 


where  J is  the  total  system  inertia  dyadic  referenced  to  point  0,  w is 
the  inertial  angular  velocity  of  a basis  fixed  in  the  spinner  (for  the 
example  considered  here),  £^  is  the  location  of  point  0 referenced  to 
the  instantaneous  cem,  and  the  vector  £ emanates  from  point  0 and 
terminates  at  a differential  mass  element. 

From  the  rotational  equivalent  of  Newton's  second  law  of  motion, 
we  have 


^d  c 

=1  (31) 

where  ^ is  a vector  of  all  external  torques  referenced  to  the  center  of 
mass  of  the  composite  system.  The  inertial  time  derivative  of  the  total 
system  angular  momentum  can  be  expressed  as 


Nj 

o 

d „c 
dt  — 

= °J  • 

W + (w  X °J  • 0))  + ° J • u 

/•  00 

X £^  m.j,  + 1 £ X £ dm 

r ° 

+ 0)  X 

I (r  X r)  dm 

J - - 

The  terms  in  Eq.  (32)  will  be  developed. 


(32) 


R.  E.  Roberson,  Dynamics  and  Control  of  Rotating  Bodies,  New  York, 
Academic  Press,  to  be  published. 
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The  integral  J Ji^  Ji  dm,  which  is  a relative  angular  momentum, 
can  be  evaluated  as  follows.  For  the  case  of  particle  masses  illustrated 
in  Figure  2,  we  have 


/o  /*  o 

r X £ dm  = r / 


or  Eq.  (33)  can  be  alternatively  expressed  as 

- o o 

/r  X r dm  = m S (r.  + 6.)  x 6. 

J-  - i -I  -I  -1 


The  derivative  of  Eq.  (34)  with  respect  to  a basis  fixed  in  the  hub 


/•  o /•  oo  oo 

"dt"  J L - J = m E (jr . + ^. ) x Aj 


The  vector  can  be  expressed  in  terms  of  deformation  of  the  particle 
masses  as  [Eq.  (4)] 


ill  « 1 

•R)  m_  . —I 


The  second  inertial  time  derivative  of  Eq.  (36)  is 


r°°  ° 

= - — E 6.  + 2wx6.+wx(wx6.) 
m„  ; L-i  - -1  - - -1 


+ « X 6 .1 


The  vector  cross-product  m^p  ^ x is  functionally  dependent  on 
products  of  and  its  derivatives;  as  such  it  can  be  neglected  in  first 
approximation  analysis. 

The  composite  inertia  dyadic  referenced  to  the  base  point  O can 
be  expressed  as 


o o 

J = I + L m.  [(r.  +_6.)  • (r.  +_6.)  E 


- (r.  +^i)  (r.  +£.)]  (38) 

where  1°  represents  the  inertia  of  the  hub  referenced  to  point  0 and  the 
summation  represents  the  contribution  of  the  flexural  subsystem  to  the 
total  system  inertia  referenced  to  point  0.  The  rotary  inertia  of  the  par- 
ticle is  assumed  negligible.  Equation  (38)  can  be  alternately  expressed 
as 


J = {x}”^  [l°  - m E (^^pTT)  (r.  + 6.)]{x}  (39) 

i 

where  the  operator  (~)  represents  a skew  symmetric  operation  on  the 
components  of  a (3  x 1)  column  matrix  to  yield 


(40) 


3Z 


The  elements  of  the  (3x1)  column  matrix  {x}  are  unit  vectors.  A first 
approximation  expression  to  Eq.  (38)  is 


°J  = }[l°  - m E r.  r,  - m S(r.  6.  + 6.  r.)]  {x} 

i i ^ 


{x  }[l  - m E(r.  + 6.  r.)]  {x} 


The  derivative  of  Eq.  (39)  with  respect  to  a platform  fixed  frame  of 
reference  is 


° I 

dt 


= J = {x}^[  - mE6.  (r.  + 6.) 

' .11  1 


- m E(r.  + 6.)  6.]  {x} 
i 


where  differentiation  has  been  carried  out  term-wise  on  only  deformation- 
dependent  terms. 

Substituting  Eqs.  (34)  and  (35)  into  Eq.  (32)  and  rewriting,  we  have, 
to  first  approximation  for  the  case  of  no  external  torque. 


, o , , o ... 

J"U+UX  J’W  + OJ'W  + X 


OO  r O T 

+ m E (_£.  + ^. ) X ^.  + w X [ m E (_r . + _6 . ) x ^.  ] = ^ 
i i 


Following  the  previously  established  pattern,  consider  the  system  to  be 
participating  in  pure  spin  in  its  nominal  equilibrium  state.  The  inertial 
angular  velocity  of  the  system  in  its  equilibrium  state  is 

n = oz 


j 
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The  inertial  angular  velocity  of  the  basis  fixed  in  the  hub  is  represented 


w = ^ + 6w 

where  the  angular  velocity  chain  rule  was  used.  A first  approximation 
expression  is  developed  by  substitution  of  Eqs.  (41),  (42),  and  (45)  into 
Eq.  (43)  and  truncation  of  terms  involving  products  of  deformation, 
deformation  rates,  variational  angular  velocities,  and  variational 
accelerations  to  yield 


x'^  [l'  6d)  + 6(»)I  f2+f2I  6co 


+ Cl  (I  - mL  6.  r.  - mS  r.  6.)  n 


+ -mE5.  r.Cl-mSr.  6.C1 

11  11 


+ m L r.  6.  + Cl  m L r^  6.]  = 0 


■^Vith  some  specialization  to  a diagonal  inertia  matrix,  the  corresponding 
scalar  equations  coordinatized  in  a platform- fixed  basis  are 


A 6u  + Cl(C-B)  6 0)  + E m.  r Me*  + 6^  = 0 

X ' y 1 y Z Z 


B 6 (1)  + Cl  (A-  C)  6 <»>  =0 

y X 

C a<j  + 2 m E6‘  n r*  - mSr*  6*  = 0 
z y y y ^ 


Assume  the  particles  are  symmetrically  located  about  the  hub  so 


r^  = -r^  = r 
y y y 
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Applying  the  transformations  [Eq.  (11)],  we  can  express  Eq.  (47)  as 


(C-B)  + mr  1 § + § ] = 0 

y / z z 

6<I»  + n (A-C)  6w  = 0 
y X 

+ ZCir^m  = 0 (49) 

The  dependent  variable  T],  which  is  associated  with  a symmetric 
mode,  is  not  present  in  the  rotational  equations.  Earlier  we  observed 
that  rotations  of  the  hub  did  not  couple  into  the  symmetric  mode  shape. 

The  excitation  of  the  symmetric  mode  shape  does  not  influence  the  hub 
rotational  motions,  nor  do  hub  rotational  motions  tend  to  excite  sym- 
metric modes  of  the  structure.  In  addition,  there  is  a decoupling 
between  in-  and  out-of- plane  motions. 

D.  SYSTEM  EQUATIONS 

The  system  equations  are  comprised  of  the  total  composite 
equations  of  motion  [ Eq.  (49)]  and  flexural  subsystem  equation  (12) 
specialized  for  the  radially  mounted  case,  repeated  here  for 
convenience: 


Out  of  Plane 


A 

B 


In  Plane 


Out  of  Plane 


A 6d>  + n (C-B)  +m  r (§  + ) = 0 

X y y z z' 


B 6(1)  + n(A-C)  6CJ  = 0 
y X 


In  Plane 


C6<1)  +2nrm5  -mF  § =0 

z y y y X 


? - 20  I + k/m  § - 26(1)  r = 0 

X y X z y 


5^  + 20  § + (k/m  - 0‘)  C - 46(*>^  ^ T = 0 

Y X.  Y z,  y 


? +(n  +k/m)?  +26w  r + 26u)  n r = 0 (50) 

z z y y 7 
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The  equations  describing  wobble  motions  decouple  as  follows: 


A 6d)  + fi  (C-B)  6oj  + m r (§  + fi  § ) = 0 

X y y z ’z 


B 6(b  + n (A-C)  60J  = 0 
y X 


+ (n  + k/m)  + 26<1)^  ^ ^^y  ~ ° 


The  additional  equations  descriptive  of  hub  motions  along  the  spin 
axis  can  be  decoupled  as  follows: 


C 6d)  + 2n  r m|  -mF  5 =0 

z y y y X 


5 - 2n  I +k/m  § - 26w  T = 0 

X y X z y 


I + 2n  5 + (k/m  - n^)  ? - 460)  n r = 0 

y X y z y 


(51) 


(52) 
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m.  STABILITY  ANALYSIS 

A stability  analysis  is  provided  for  the  simplified  spacecraft  system. 
Emphasis  is  on  the  development  of  literal  stability  criteria  of  value  in  pre- 
liminary design  phases.  Admittedly,  some  risk  is  associated  with  extra- 
polation of  results  for  simplified  dynamical  systems  to  more  complicated 
examples.  The  stability  criteria  can  be  of  value  in  preliminary  design 
phases  and  are  not  intended  to  supplant  more  detailed  followup  studies  to 
verify  ultimate  system  stability. 

The  results  presented  in  this  section  parallel  closely  the  work  of 
P.  Likins  and  F.  Barbera.  ^ Their  work  addressed  the  problem  of  develop- 
ing literal  attitude  stability  criteria  for  idealized  spinning  spacecraft.  The 
stability  criteria  essentially  represented  a duplication  of  results  from  sev- 
eral prior  studies  involving  considerably  more  sophisticated  models  as  a 
point  of  departure.  The  essential  resvilts  of  that  effort  are  described. 

A.  STABILITY  ANALYSIS  OF  SIMPLE  PARTICLE  MODEL 

The  equations  descriptive  of  the  wobble  motion  for  the  simple  particle 

model  of  Figure  1,  provided  in  Section  11  [Eq.  (5l]  , are  repeated  here  for 

* 

convenience: 

A 6(b  -(B-C)n6w  +mr  (F  ) = 0 

x'  y y ^z  ^z 

B + (A-C)  =0 

y X 

’e'  + 2Q  Of  + t + 26cb  F +2  6(o  OF  =0  (53) 

bz  -z  ^z  X y y y 

3{c  ^ 

An  additional  damping  term  has  been  included  in  the  third  equation  {Q  is  the 
damping  ratio),  a represents  the  loaded  frequency  for  the  radially  mounted 
subsystem. 
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A,  B,  and  C represent  total  system  inertias  about  x,  y,  and  z respectively; 


f is  the  difference  between  the  z axis  deformation  of  the  particles;  and 
• z 


6oi3^  and  are  the  angular  velocity  variational  coordinates  from  nominal 


spin.  The  stability  of  that  portion  of  the  system  descriptive  of  wobble  mo- 
tions can  be  obtained  by  examination  of  the  characteristic  equation  using 
Routh-Hurwitz  criteria,  which  yield  both  necessary  and  sufficient  conditions. 

The  characteristic  equation  associated  with  Eqs.  (53)  can  be  expressed 


as 


+ [ 2C  cr  + 2C  o K'] 


+ [ + K'  (a^  - £2^] 


+ [ 2£  a + 2C  a K'  a^]  S 


+ [ + K'  (a^  - n ^)]  =0 


(54) 


where 


2m  r 

^ 


A' 


^2  (C^  - A') 

B' 


2 _ (C"  - B') 


& = 


A' 
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I 


and  A^  B',  and  represent  core  inertias.  The  characteristic  equation  is 
of  the  following  form: 


+ P3  + P2  S + Pq  = 0 


(55) 


One  form  of  the  Routh-Hurwitz  criteria  yields  the  following  necessary  and 
sufficient  conditions  for  asymptotic  stability: 


Criterion  1 p^  > 0 


Criterion  2 p<  > 0 


2 2 

Criterion  i p^  P2  P3  " Pj  " Pq  P3  ^ ® 


Criterion  4 p^  > 0 


The  results  obtained  by  applying  these  criteria  for  the  system  under  consid- 
eration are  summarized  in  Table  3. 


Criterion  4 can  alternatively  be  expressed  in  terms  of  hub  inertias  B C ^ 


as 


,.Z  2m 

(1)  > 1. 

2mr^  + C"-B' 

y 


(56) 


The  loaded  frequency  a for  the  case  of  the  radially  mounted  particle,  given 
in  Table  1,  is  repeated  here  as 


Jx>aded 


k/m  + S2 
Natural 


(57) 


Frequency  Frequency 


Preload 
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Si 

Table  3.  Necessary  Conditions  for  Stability  of  Spacecraft  Motions 


Criterion 

Condition 

Comments 

1 

£ > 0 

Positive  damping 

2,3 

C>  A 

OB 

Spin  about  cixis  of 
maximum  moment  of 
inertia  (system) 

4 

2 2m 

(-)  > ^ 

'fl'  ^ C - B 

Bounding  criteria  on 
loaded  natural  fre- 
quency in  terms  of 
spin  frequency  and 
system  inertia 
properties 

^The  criteria  are  both  necessary  and  sufficient  conditions  for  asymptotic 
stability  for  equations  associated  with  the  wobble  motions  of  the  space- 
craft. The  condiUcns  are  only  necessary  conditions  for  asymptotic 
stability  of  equations  of  spacecraft. 

and  is  augmented  positively  by  the  effects  of  preload.  If  the  hub  spin  inertia 
exceeds  the  transverse  inertia,  then  the  inequality  [Eq.  (56)]  is  trivially 
satisfied  and  furthermore  is  independent  of  the  amount  of  flexibility  or  the 
magnitude  of  spinfl. 

The  dynamical  model  and  the  associated  stability  criteria  developed  for 
the  wobble  motion  of  the  simplified  spacecraft  system  are  sufficiently  general 
to  permit  examination  of  alternative  ways  of  attaching  the  flexural  subsystem. 
We  observed  in  Section  II  that  the  anti  cantilever -mounted  substructure  ex- 
perienced a reduction  in  its  natural  frequency  due  to  the  preload  effect,  whereas 
the  intermediate  case  did  not  experience  a pulling  of  its  natural  frequency  due 
to  preload  effects.  If  we  assvime  that  the  hub  inertias 


C'  = B* 
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then  the  anticantilever -mounted  subsystem  would  lead  to  more  stringent 
criteria  on  the  unloaded  structural  subsystem  to  ensure  that  system  wobble 
motions  are  stable.  We  can  easily  see  this  by  substituting  the  expression  for 
the  loaded  frequency  of  the  anticantilever  subsystem 

cr^  = k/m  - ( 58) 


into  the  inequality  [Eq.  (56)]  and  observing  that  stability  is  ensured  only  if 

the  unloaded  structural  frequency  squared  is  at  least  twice  as  great  as  the 

2 

spin  frequency  squared  (k/m>  2S2  ).  The  orthogonally  mounted  subsystem 

requires  that  the  unloaded  structural  frequency  be  greater  than  the  spin 
1/2 

frequency  [(k/m)  > ] . 

B.  COMPARISON  WITH  PREVIOUS  STUDIES 

The  literal  stability  criteria  presented  have  been  developed  for  a simple 
particle  construct  that  is  supposedly  representative  of  a more  general  flex- 
ural subsystem.  The  fidelity  of  such  criteria  in  predicting  stability  require- 
ments for  more  complex  spacecraft  might  well  be  questioned.  Likins  and 
Barbera  provided  comparisons  with  stability  analysis  results  published  in  the 

literature  for  a specific  class  of  spacecraft  exhibiting  flexural  booms  movmted 

7 - 10  . 

radially  outward  from  a rigid  core.  With  proper  interpretation  of  results. 


*7 
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<1 


the  literal  stability  criterion  developed  for  the  simple  parti  :le  idealization 
and  the  stability  criterion  developed  for  models  which  were  more  complex 
in  idealization  give  surprisingly  close  agreement.  The  details  of  these  com- 
parisons are  discussed. 

•7 

Rakowski  and  Renard  studied  the  nutational  behavior  of  a torque -free 
spinning  satellite  composed  of  a rigid  hub  and  flexible  appendages  in  an  equa- 
torial plane.  The  vehicle  studied  is  referred  to  in  the  literature  as  a crossed 
dipole  configuration.  The  flexible  appendages  are  two  pairs  of  uniform 
booms  located  along  the  trauisverse  axes  of  principal  inertia  of  the  hub.  Two 
separate  classes  of  antenna  flexure  considered  were  equatorial  plane  and 
meridian  plane  motions.  The  following  limiting  assumptions  were  made: 
small  displacement  vibrations,  nearly  constant  spin  speed  of  the  satellite, 
and  small  products  of  the  transverse  angular  rates  as  compared  to  the  square 
of  the  boom  natural  frequency.  These  assumptions  validate  the  use  of  pre- 
viously determined  natural  frequency  and  mode  shape  data  that  reflect  the 
effects  of  spin.  The  fundamental  mode  of  vibration  was  assumed  to  predomi- 
nate. The  justification  for  the  truncation  of  higher  modes  was  that  the  rigid 
body  forcing  occurred  at  a frequency  closest  to  the  loaded  lowest  natural  fre- 
quency of  the  booms.  The  response  in  the  higher  modes  to  that  excitation  was 
deemed  small.  The  stability  of  this  configuration  for  a variety  of  system 
parameter  variations  was  assessed  by  numerical  integration  of  the  full  non- 
linear equations  of  motion.  Stability  boundaries  were  developed  in  plots  of 
normalized  inertia  properties  for  different  values  of  a nondimensional  param- 
eter that  reflects  a ratio  of  centrifugal-to-elastic  forces.  These  parameters 
are  defined  in  terms  of  variables  used  herein.  Let  the  normalized  inertia 
parameters  be  defined  as  I 


where  Ig  is  the  inertia  of  a unifornaly  distributed  beam  about  the  core. 


In 


the  simple  particle  model,  the  analogy  would  be 


(61) 


The  nondimensional  ratio  \ is  defined  as 


El/pf 


(62) 


A necessary  condition  for  stability  for  a particle  mass  approximation 
to  a crossed  dipole  configuration  is 


(^) 


2m.  d^ 
1 1 


2m.  d . + - I. 

11  1 


(63) 


.th 


where  o).  is  the  loaded  frequency  associated  with  the  i axis,  d.  is  the  nomi- 
1 ^2^  ^ 
nal  displacement  of  the  i particle  from  the  axis  of  spin,  and  I is  the  inertia 

.tb 


of  the  hub  about  the  i"*'  axis.  The  necessary  conditions  for  stability  for  the 
crossed  dipole  configuration  can  be  broken  down  into  two  separate  criteria 
that  are  nearly  identical  with  the  criteria  developed  for  the  case  of  a rigid 
hub  and  a single  pair  of  radially  mounted  booms. 


Developed  in  Ref.  5 and  summarized  in  the  Appendix 


r. 

/■' 

r 
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The  only  other  identification  we  must  make  to  complete  the 
correspondence  of  the  literal  stability  criteria  results  with  the  work  of 
Rakowski  and  Renard  is  the  loaded  frequency.  The  loaded  frequency  of  a 
uniform  boom  is^’  ® 


Z aZ  , j 
(0  = a + 1 . 193  fi 


(64) 


where  5 is  the  lowest  unloaded  natural  frequency  of  a uniform  boom  and 
11 

IS  given  as 


0=  (3.  515) 


El 


1/2 


(65) 


By  comparison,  the  loaded  frequency  of  a simple  particle  analogy  is 

Z aZ  , iLL\ 

o)  = a + S2  (66) 

1 

I Thus,  despite  a significant  difference  in  models,  the  loaded  frequencies 

are  nearly  the  same. 

The  simple  particle  stability  criteria,  with  minor  interpretation,  can 
be  reexpressed  in  terms  of  the  nondimensionalized  parameters  used  in 
Rakowski  and  Renard' s work.  The  original  necessary  condition  for  stability 
can  be  reexpressed  as 

h 

1 

I 

I 


I 


^ Z 2 Z 

Likins  and  Barbara  used  a loaded  frequency  of  = a*"  + 1 . Z , which  is 

generated  for  a massless  uniform  beam  with  particle  tip  mass. 

^W.  C.  Hurty  and  M.  F.  Rubinstein,  Dynamics  of  Structures,  New  Jersey, 
Prentice -Hall,  1964. 
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(67) 


S2  pi 


+ 1 - (i/Kp) 


or 


— r 

+ l.  193  > 


fr  + 1 - (i/Kp) 


(68) 


Rakowski  and  Renard's  results  are  reproduced  in  Figure  7.  The  results  of 
Barbera  and  Likins,  which  amount  to  a literal  stability  criterion  developed 
for  a simple  particle  analogy,  compare  very  favorably  after  some  minor 
reinterpretation  outlined  above. 

Additional  comparisons  have  been  provided  in  Ref.  5 with  other  anal- 
8 10 

yses.  ’ These  comparisons  are  not  reiterated  here.  All  these  compari- 
sons substantiate  the  fact  that  the  results  developed  for  the  particle  analogy, 
with  some  minor  reinterpretation,  can  serve  as  a useful  starting  point  in  pre- 
liminary design  work.  As  in  most  studies,  there  are  certain  limiting-case 
situations  for  which  the  criteria  might  be  misleading:  for  example,  high  spin- 
to-stiffness  ratios.  At  the  limiting  case  of  high  spin -to -stiffness  ratios,  we 
would  expect  that,  as  the  necessary  condition  for  stability,  the  hub  spin  inertia 

must  exceed  a transverse  inertia  of  the  hub.  And  yet,  if  we  use  the  approx- 

2 ..2  2 

imation  for  the  loaded  frequency  u)  = {o  +1.20),  then  the  requirement 
(C*  > B')  does  not  emerge  from  the  stability  ineqxiality.  At  high  spin-to- 
stiffness  ratios,  the  loaded  frequency  is  close  to  the  loaded  frequency  pre- 

2 A 2 2 

dieted  for  the  simple  particle  anadogy  (w  = a +f2  ).  The  ratio  of  loaded 

frequencies  to  spin  frequency  for  a uniform  massless  beam  with  a particle 

5 

tip  mass  is  analytically  given  as 


* 


This  result  can  be  verified  in  Ref. 


9. 
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k 12 


\ 

\ 

L 


a (QUASI -RIGID  BODY) 
b X = 10 
c X = 100 

O RESULTS  OF  LIKINS  AND  BARBERA 


A 

r = g 

c’ 

^ A C 
P B’  . 

A 

El/pf 


(69) 


“ ~ |tanh  [(3)^^^  (J2/a)l  {n/5)l| 

and  can  be  verified  to  approach  a value  of  1 as  S2/d  approaches  infinity.  Thus, 
the  discrepancy  arises  in  approximating  the  loaded  frequency  of  the  uniform 
massless  beam  with  particle  tip  mass  as 


2 A 2 


+ 1.2 


(70) 


as  opposed  to  its  true  value  at  high  spin -to -stiffness  ratios,  which  in  the 

....  . .(i>  . 12  . 

limiting  case  is  one  (?r  -•  1 as  — -♦  ®). 

a 

C.  STABILITY  CRITERIA  EXTENSIONS 

To  date,  literal  stability  criteria  have  been  developed  for  simplified 
particle  mass  idealizations  of  a spacecraft.  Literal  stability  criteria  have 
been  developed  for  more  general  structures  configured  such  that  the  append- 
age particles  in  their  undeformed  state  lie  in  a plane  perpendicular  to  the 

5 

principal  axis  of  spin  and  passing  through  the  system  center  of  mass. 

These  developments  are  algebraically  cumbersome  and  do  not  permit  the 
clear  interpretations  provided  in  the  simpler  examples.  As  a result,  only 
the  basic  outline  of  the  approach  taken  is  provided,  along  with  the  major 

5 

results  of  the  study.  The  reader  is  referred  to  the  original  study  for  a 
more  complete  explanation. 

The  equations  of  motion  for  the  appendage  are  developed  for  the  i^^ 
particle  of  the  appendages  idealized  as  described  above.  The  equation  devel- 
opment is  nearly  identical  to  the  development  in  Section  II.  An  equally 
acceptable  starting  point  might  be  simplification  of  the  more  general  structure 
of  equations  provided  in  Ref.  1.  Further  assuming  that  the  contact  forces 
resolved  in  the  transverse  plane  are  independent  of  z axis  deformations  or 
their  rates  or  of  contact  forces  resolved  along  the  z axis  deformations  or 
their  rates,  we  can  verify  the  wobble  axis  motions  do  decouple  from  in-plane 
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motions,  as  observed  previously  for  the  limiting  physical  example.  The 
resulting  equations  contain  the  dependent  variables  ^ axis  deforma- 

tion of  the  i*^^  particle,  and  the  two  transverse  variational  angular  velocities 

CO  and  oo  . Collectively,  these  generic  equations  for  the  appendage  subsystem 
X y 

are  given  in  Eq.  (7  1 ) in  a matrix  format  with  the  vector  q identified  as  an  N 
vector  of  z axis  displacements.  The  appendage  subsystem  equations  reduce 
to 

M q + Kq  = M T f (<I)  - w ) - MF  I (w  + flo)  ) (71) 

X y X y X y 

where  the  contact  forces  have  been  replaced  by  their  matrix  counterparts, 
and  where 


M = 


m 


m. 


1 M, 


-m.m 
1 n 


-m,  m 
1 n 


m 

I 

n Mn 


r"  T- 1 


r = 

X 


.N 


I = 


(N  X 1) 


1 
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w 


The  wobble  motion  portion  of  the  system  equations  reduces  to 

j—  T 

A6i  - n6u)  (B-C)  + £2^  (M  r f q + (M  r f ) q=0 

X y y ' 

*1* 

B6i  +n6u)  (A-C)  -£2^  (M  r f q - (M  r f ) q=0 

y X ^ 


(72) 


The  column  matrix  q will  be  replaced  by  «T1,  where  $ is  a matrix  of  eigen- 
vectors obtained  from  the  eigenvalue  problem  associated  with  the  homo- 
geneous matrix  equation 


Mq  + Kq  = 0 (73) 

With  appropriate  normalization  of  f such  that 

j'^MS  = E (NxN)  (74) 

and  with  the  addition  of  modal  damping,  the  appendage  equations  reduce  to 


T]  + 2;  «T1  + w T)  = 6 (6i  - nfitt)  ) 


6 (6i  + 06(1)  ) 

y X y 


(75) 


*At  the  risk  of  duplicating  nomenclature,  the  variable  T)  was  selected  to 
represent  a modal  deformation.  The  same  variable  was  used  earlier  in 
the  text,  the  definition  of  which  appears  in  Eq.  (il). 
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where 

6 = M r i (N  X 1) 

X X 

6 = M r i (N  X 1) 

y y 

and  the  equations  descriptive  of  wobble  motion  reduce  to 
A6i^  - n6(Uy(B-C)  6y  T1  + 6y  T1  = 0 

B6i  + neiu  (A-C)  - = ® 

y X AX 

These  equations  are  subsequently  specialized  by  truncating  to  one  modal 

coordinate.  The  (N  x 1)  vectors  6 and  6 reduce  to  scalars  designated  as 

X y 

J^^and  6yj,  With  the  corresponding  reduction  in  dimensionality,  literal 
stability  criteria  can  be  developed  by  Routh-Hurwitz  analysis.  The  neces- 
sary and  sufficient  conditions  for  asymptotic  stability  of  the  wobble  motion 
are  given  in  Table  4,  based  on  the  indicated  restrictions. 

Table  4.  Necessary  and  Sufficient  Conditions  for  Stability 
of  Wobble  Motion  for  Single  Mode  Case 


C>A  1 

>•  Major  axis  spin  requirement 

OB  J 

* *y, 

[irj  ^ (C-A)  (C-B) 
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In  the  previous  analysis,  one  restriction  was  that  all  modes  were 
truncated  with  the  exception  of  a single  modal  frequency.  Can  satisfaction 
of  the  inequality  in  Table  4 for  each  mode  individually  ensure  stability  for 
the  system  in  which  all  modes  are  present?  Necessary  conditions  can  be 
developed  using  Routh-Hurwitz  analysis  for  the  multimode  case.  The  devel- 
opment can  be  algebraically  encumbered  if  we  attempt  to  develop  all  the 
necessary  and  sufficient  conditions  for  stability.  However,  by  concentrating 
solely  on  a single  inequality  of  the  Routh-Hurwitz  criteria,  necessary  condi- 
tions can  be  identified.  In  particular,  the  requirement  that  the  constant  ter 
from  the  characteristic  equation  be  greater  than  zero  [Eq.  (55)]  yields  the 
criteria  indicated  in  Table  5 for  the  cases  indicated. 

Table  5.  Necessary  Conditions  for  Stability  of  Wobble 
Motion  for  Multimode  Case 


r = 0 

X 


T 2 

6 (<*>  ) 6 


r = 0 
y 


r , r ^ 0 

X y 


T T 

6 6 -6  6 =0 
X y y X 


T 2 ‘ 


, (C-A)  6j  6„  + (C-B)  6^ 

(C-A)  (C-B) 


The  subject  of  enhancing  stability  of  a spinning  flexible  body  through  active 
control  was  briefly  discussed  in  Ref.  5.  The  results  achieved  to  date  are 
summarized  in  the  Appendix.  This  area  in  particular  deserves  considerable 
attention  and  will  be  the  subject  of  much  continuing  investigation. 
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APPENDIX 


SUMMARY  OF  STABILITY  CRITERIA 


The  following  criteria  have  been  formally  established  as  necessary 
**  . 


conditions  for  asymptotic  stability  of  spin. 

C>  A AND  C>  B IN  ALL  MODELS 


MODEL 


CRITERION 


2 j2 

'n'  (C-B) 


2 , 2 
'n'  (c-A) 


2 o j2 

0).  Z md 
kl'  (C-B) 


[ 


Reproduced  in  its  entirety  from  Ref.  5 by  permission  of  the  authors. 
In  some  cases  sufficiency  has  also  been  proven  formally. 
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MODEL 


CRITERION 


z 


2 


o « 2 
2 m jf 

(C-A) 


0 


2m  j d 


2m  j d 


+ C'  - B' 


2m2i 


2 m2i 


+ C"  - A" 


w~  loaded  frequency  associated  with  m^ 
a~  loaded  frequency  associated  with  m^ 


w ^ (C-A)  m d^  + (C-B)  m I ^ 

'd'  (C-A)  (C-B) 


54 


MODEL 


CRITERION 


r,  = r,  = o 


r = r = o 

1 


r,=o 


- .T  r Zl-l  . 

.2  6 l(<)  J 6 

(1)  >-J£— : 

'n'  (c-B) 


2 sj  [,.2]-'  e 

© > ^-(C-A) 


Mode 


co^  ^ (C-A)  + (C-B) 


<?f)  > 


(C-A)  (C-B) 


Mode  (6  6^-6  dj)  = 0 

X y y X 


(C-A)  6^  6 + (C-B)  [</]'^  6 


(C-A)  (C-B) 


<1. 


MODEL 


CRITERION 


MODEL 


CRITEPION 


k = 00 

y 


m 


(j)  2 

(— ) > 

'fi  ' ^ (C-A) 


4m 


u_  2 m (C-A)  + (C-B) 


(jf)  > 


r cm^r^r^ 
y 


(a)  2 -ifc 

LC(Tf) 


where 


D = (C-B) 


CmT' 


(C-A)  - 


CO  2 _ 

C (ji)  + 4m  r 


/V 

z 


First  Mode 


C*t>AandC+^>B 


(0 . 2 
(cT)  > 


xl 

17“ 


SYMBOLS 


Operator 

-4  ( ) = (") 

dt 


■ar 


‘^  ( ) = (°) 


(■") 


Scalars 
A,  B,  C 

A',  B',  C' 


o 


a 


second  inertial  time  derivative  of  quantity  in  parentheses 


time  derivative  of  quantity  in  parentheses  with  respect 
to  a rotating  frame  of  reference  P 

skew  symmetric  operator  that  transforms  a (3  X 1)  matrix 
into  a (3  X 3)  matrix;  Eq.  (16) 


transpose  operation  on  a matrix 


principal  inertias  of  system  in  equilibrium  with  respect 
to  point  0 referenced  to  the  x,  y,  and  z axes  respectively 

principal  inertias  of  rigid  portion  with  respect  to  point  0 
referenced  to  the  x,  y,  and  z axes  respectively 

unloaded  frequency  of  oscillation  reflecting  elastic  restor- 
ing force  effects  only 

mass  of  particle  i 

total  mass  of  spacecraft 

mass  per  unit  length  of  uniform  boom 

distance  of  particle  in  its  equilibrium  orientation  from 
point  0 (rl  = y) 

loaded  frequency  for  radially  mounted  subsystem 
(k/m  +n2)l/2 

lowest  unloaded  natural  frequency  of  uniform  boom 


SYMBOLS  (Continued) 


Vectors 


S 

j -J 

F'^ 


H 


c 


r. 
— 1 


6. 

— 1 


JL 

1 

— 

6cJ 

n 


portion  of  contact  forces  applied  to  particle  i due  to 
effects  of  elasticity 

vector  sum  of  i^ll  contact  forces  applied  to  particle  i (j  is  an 
arbitrary  summation  index) 

portion  of  contact  forces  applied  to  particle  i due  to  effects 
of  spin;  also  known  as  preload  effect 

total  angular  momentum  of  composite  system  referenced 
to  composite  center  of  mass 

nominal  equilibrium  position  of  particle  i with  respect 
to  point  0 

instantaneous  position  of  particle  i with  respect  to 
Newtonian  frame  of  reference 

small  deformation  of  particle  i from  its  equilibrium 
position 

vector  sum  of  particle  deformations  from  equilibrium 

vector  emanating  from  ccm  and  terminating  at  point  0 

vector  difference  of  particle  deformations  from 
equilibrium 

inertial  angular  velocity  of  P frame  of  reference 

variational  angular  velocity  from  pure  spin;  Eq.  (9) 

nominal  inertial  angular  velocity  of  (x,  y,  z}  when  system 
is  in  equilibrium  orientation;  assumed  to  be  a state  of 
constant  spin  (Q  = Oz) 


Dyadics 


E 


unit  dyadic  (E  = + yy  + zz) 


SYMBOLS  (Continued) 


